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1, INTRODUCTION 

Most modern high performance radars use travelling-wave tube amplifiers to obtain coherent 
transmission [1]. As pointed out previously, these tubes work more efficiently under constant amplitude 
conditions. Moreover, good amplitude modulation (AM) is difficult (and very expensive) to achieve with 
these devices. Therefore, only purely phase modulated pulse trains have been considered so far. However, 
with the emergence of solid state microwave sources, efforts are being made to replace the relatively large 
and expensive vacuum devices by low power solid state elements and the waveguide elements by planar 
circuits. With these new components the size and cost are reduced and a number of commercial applications 
become feasible [2]. 

The contribution of solid-state devices has also been significant in areas where performance was 
previously inadequate. For example, it is much easier to use any form of modulation with solid state 
components [3]. Although AM is not an efficient method to provide the large time-band width required for 
good radar performance, [4,5] it does provide an excellent means of improving the resolution capability. 
Consequently, this chapter treats the problem of finding ER amplitude and phase modulated (a.m and p.m) or 
multilevel sequence. 
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2. CHOICE OF STARTING SEQUENCE 

The choice of the starting sequence is an important process to determine the best sequence of 
random variables, describing an outcome following the deterministic pattern and probability distributions. 
These constructs are extremely useful in probability theory. 


2.1. Choice at Random 

Random sequence with optional properties that use property such as integer or length value are able 
to specify in any methods and any order to get the possibilities sequence intended for serving as a starting 
sequence. 

The design of a Huffman sequence in general requires the choice of the code length (N+1), the circle 
radius X and the zero pattern for which the magnitude sequence la(n)l is most uniformly distributed. 
However, it remains as an unsolved problem. The direct evaluation of all 2N possible root pattern for a given 
radius X and N 1s not feasible if N is larger (N>20). It remains invariant even after account is taken of the 
zero patterns formed from others by rotation in the Z-plane through an angle @ or by other transformations. A 
general trial of the error procedure was done to chooe a root pattern randomly, to compute the sequence for a 
succession of values for radius X [6,7]. 

This method focuses more into this problem including the stationary phase principle applicable for 
an arbitrary length sequence. Zero pattern is one of the random solution on this method. In fact, Huffman, 
1962 suggested that in order to maximize the ER, the roots should be chosen in a random fashion with 
approximately half the zeros on each circle. He further concluded that the optimum ER should be 
proportional to (N+1)-1/2. Generally known that the correlation properties are used to judge the random 
sequence whereby the sequence is uncorrelated with itself, i.e. random, if its ACF has uniformly low side 
lobes. Although various methods as described above is acceptable for designing Huffman codes with 
moderately large length (N=100), the major difficulty arises for longer sequences because of the difficulty of 
the computational nature due to factorizing polynomials of degrees which are larger than 100. The following 
sections describe some useful methods for designing the multilevel sequences [8]. 


2.2. Designing Multilevel Sequences using m integer and N length 

To generate the multilevel sequences the length and some integer is considered. The impulse- 
equivalent sequence is a finite sequence of complex numbers whose autocorrelation is zero +N [9]. The zero 
patterns in Figure 2.3 had been discussed in Chapter 2. The zeros of +N lies on 2 circles is an impulse 
equivalent which is a complex number of finite sequences. 
a) Select m with any integer (m=1, 2,3...) and any length of ‘N’. 
b) X-X-1 should be an integer coefficient and ~ is an odd integer and therefore choose X to be the positive root of X- 


X-1. 
a’ —bX —c=0 (1) 
ya btvb7-4ac 


2a (2) 
Equation | is a quadratic equation for: Using ‘quadratics coefficient’, applying this formula (Equation 3). 


X* —mxX —1=0 (3) 


a) Generated Ra Me 
b) Using 4 conditions (Equation 4,5,6 and 7) for generated integer values following length N: 
For length I until 3 














H (1) =1, H(2) =2m, H(3) =2M’ (4) 
For length 4 until N+1 
The centre element of [| = ©) +] 

H(k)=mH(4)+ H42-)-Ad*(xX £- XxX *) (5) 
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For length N is K > ~ + 1 
H(L—kN) =-[H(Kn+)](—-)*” (6) 


If else of 3 conditions above 


A(k) =mA Ck —-1)+ Hk — 2) (7) 


Table | shows is the set of sequences with 13-length of sequences and m=1, 2, 3, and 4 designing by 
multilevel sequences generated by m integer and N length produces m(x) sequences. 


Table 1. Sequences Generated by Integer Values of m 


L m Sequence S.L.E E.R P.S.E 
13 1 [12246 10 -1.89 10 -6 4 -2 2 -1] 1.1723 3.209) 0.6880 
2 [1 4 8 20 48 116 82.01 116 -48 20 -8 4 -1] 0.9723 2.9138 0.4856 
3 [1 6 18 60 198 654 862 654 -198 60 -18 6 -1] 1.0911 2.2674 0.6692 
4 [1 8 32 136 576 2440 4558 2440 -576 136 -328 -1] 0.9703 1.6070 0.6663 


2.3. Zeros Pattern 

The zeros pattern of the z-transform in generating the sequence indicates a possible position of roots 
at the radius which the circle lies. As an alternative of having a good ACF properties, the zeros pattern can be 
edited by adding new zeros with the same angle in any length of the existing zero plot. Therefore, this 
procedure can be adopted to test the stability and for analyzing the systems by inserting the multiple methods. 

Modulating pulse is achieved by transmitting the pulse that correlated the received signal. Radar, 
sonar and echography applications basically used in signal processing method by pulse compression is to 
increase the range solution and to reduce the noise [10]. 

This method is a procedure for obtaining the set of sequences by zeros plot movement using 
fdatools. The multilevel sequence also used as the initial sequence will be represented as zeros of the 
polynomial. This method is aiming to estimate the right place of the new coordinated accordance with the 
characteristic from the initial sequence. Figure 1 shows the flowchart for this procedure. 

The probability of choosing the zero elements is different from the choosing a nonzero element. The 
sufficiently large field size whether the sequence is singular or not is determined by probability one by the 
zero pattern of the sequences, 1.e. the zeros are located in zeros plot. To create the purely real multilevel 
sequence that have a poor ER or large ACF of side lobes, the plotting should lie either of two circles which 
one circle at radius X while the other radius at X-1. In order to design the multilevel sequences, there are 
three types of problems that should be faced: 

a) Firstly, choose the number of pulses in the sequence where the number of roots following by the 
algorithm varies the positions. 

b) Secondly, the radius of one or the other circle must be chosen in the complex plane on which the root 
lies and the maximum values of ER are determined. This is to ensure that it increases exponentially with 
the number of root positions which have to be investigated. 

c) Lastly, the roots must be decided on which circle that each root should lie. 

The zeros plot on the z-plane are linear system. Matlab software used to generate the program 
function to view the pole-zero plot for this filter, 1t can be use as an example: 

>>zplane (zer,pol) 

To access to an additional tools, fvtool can be used. Firstly, the poles and zeros should be converted 
to the transfer function form and then the fvtool can be used. 

>>[b,a]=zp2tf (zer, pol, 1), where ‘1’ is a gain 

Designing Multilevel Sequences using m integer and N length are produced to obtain m(x) 
sequences. Table 1 shows the examples of some integer of Huffman sequences. It is generated by m integer 
and length of thirteen, the contents of that table also include the SLE, ER and PSL. Referring from that data, 
the best value of that parameter was found. Figure 2 shows zeros plot pattern movement for length N=13 


with the integer of m equivalent to one. The length of N should be an odd integer to generate N lengths 
sequences are represented as zeros of the polynomials [11]. 
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Select m, where ‘m7’ is any integer 
(Multilevel code/sequence that have 

good Autocorrelation Function’ using 
Some Integer Huffman Sequences’) 









Choose X to be the positive root. This is for 
expressible as a sum of powers of (X-X"') 












Generated (x(2) — x-(2) for various of 
length N 


Find the Huffman sequence generated by 
integer values of ‘m’ 


No 


H(K) = mH(k —1) + H(k - 2) 





If H(k) is the 
set of N-length 
sequence? 











No 
Y (after convert from new sea) 


e 
Normalized the sequence 


Normalized, evaluate the sequences to get the SLE and 
the ER 
Editing/adding zeros poles 
Convert coordinates to polynomial sequence 


Figure 1. Flow chart of code plot coordinates with adding/editing zeros 


Original sequences generated by this methods produced H(x)= [1 2 2 4 6 10 -1.8885 10 -6 4 -2 2 -1], 
while SLE and ER is 1.1723 and 3.2557. 
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Figure 2. Zeros pattern for length-13 and integer m=1 


Using fdatools (Matlab software) to import filter from workspace to the numerator and denominator 
to enter that particular value at Pole Zero Editor. As a theory, when the z-plane of the zeros lies in the two 
balancing circle, the envelope of the pulse train tends to be relatively symmetrical in the mid-pulse [12]. All 
possible positions of the roots in the Z-domain are generated and each position which the radius circle lies 
has been optimized. The result of real and complex from the original sequences are presented in Table 2. 
That table shows the rectangular and polar coordinates from original and editing pole defined to distinguish 
both results which is the best result discovered. Figure 3 shows the complex envelopes of the pulse train from 
the original and new editing from the thirteen zero patterns. 


Table 2. Rectangular Coordinates and Polar Coordinates from Original and Editing/Adding Pole 


Rectangular Coordinates Rectangular Coordinates 


Real Imaginary Real Imaginary 

Original pole Editing/adding pole 
-2.068 0 -1.0953 0 

-1.0854 1.2063 -0.7553 0.8175 
-1.0854 - 1.2063 -0.7553 -0.8175 
0.7765 1.2366 0.5789 0.9090 
0.7765 - 1.2366 0.5789 -0.9090 
0.0425 0.9428 0.04261 0.9875 
0.0425 - 0.9428 0.04261 -0.9875 
-0.3955 0.5588 -0.5198 0.7652 
-0.3955 - 0.5588 -0.5198 -0.7652 
0.4483 0.4642 0.6312 0.6605 
0.4483 - 0.4642 0.6312 -0.6605 
0.4953 0 0.8928 0 
Polar Coordinates Polar Coordinates 
Magnitude Angle Magnitude Angle (rad) 
(rad) 

Original pole Editing/adding pole 
2.068 3.1416 1.0953 3.1416 
16227 2.3035 1.113 2.3166 
16227 -2.3035 1.113 -2.3166 
1.4602 1.0101 1.0777 1.0038 
1.4602 -1.0101 1.0777 -1.0038 
0.9438 1.5257 0.9884 5277 
0.9438 -1.5257 0.9884 -1.5277 
0.6847 2.1867 0.9250 2.1675 
0.6847 -2.1867 0.9250 -2.1675 
0.6453 0.8029 0.9136 0.8081 
0.6453 -0.8029 0.9136 -0.8081 
0.4953 0 0.8928 0 
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Figure 3. Complex Envelope of the Pulse Trains from the Original and editing Zeros pattern 


3. RESULTS AND DISCUSSION 
3.1. Effect of Optimum Pulse Position 

In the original polynomials where the function of rectangular coordinates contains the zero-pole- 
gain such as poly (rec_coordinates) function. Figure 3 shows the movement of the root coordinates which 
tends to lie closer in the unit circle. That approached was done by implementing zeros plot (adding and 
editing) which has been proven effective in creating new illustrations. All of these zeros will be decided on 
their own about where the root in that circle should lie. Reference from two different sequences is a 
computational procedure which gives accurate result even when the SLE is low while the ER is mostly high. 
From the adding and editing method, real and imaginary axis from rectangular coordinates Table 1 illustrated 
as ‘p’ for evaluation by Matlab using the poly(p) function where p is a vector to return a row vector whose 
elements are the coefficients of the polynomial. With that, it produced the new set of sequence H(b) = 
[1.0000 0.2473 0.9478 0.4864 1.3280 0.3198 0.0088 0.5767 -1.1750 0.4685 -0.9392 0.3623 -0.9817]. The 
SLE of H (b) sequence is 1.4635, while the ER is 4.5202. Figure 4 and Figure 5 shows the autocorrelation 
function values of both sequences. Both of this graph have a positive result individually where the Hx 
sequences have the low SLE while the Hb have the best PSL and the details of the value is in Table 3 and 
Table 4. The original sequences H(x)= [1 2 2 4 6 10 -1.8885 10 -6 4 -2 2 -1] 
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Figure 4. Autocorrelation Function of the Original Sequences (Hx) 


Table 3. Correlation Values Obtained within input Sequences of Length 13 (Hx) 


T9 hx "Y hx "2 hx 13 hx any Ts hx YG hx "7 hx Yg hx Yo hx "10 hx 11 hx 


1 -0.1160 0.6880 —-0.0464 0.2949 = -0.0232 0.0983 — -0.0003 0 -0.0024 -0.0042 -0.0049 


The generated sequence H(b) = [1.0000 0.2473 0.9478 0.4864 1.3280 0.3198 0.0088 0.5767 -1.1750 
0.4685 -0.9392 0.3623 -0.9817] 
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Figure 5. Autocorrelation Function of the Generated Sequences (Hb) 


Table 4. Correlation Values Obtained within input Sequences of Length 13 (Hb) 
"Ocnb) "1 hb) 2H) 3 (ab) Tanp) TS (nb) 6a) "7 ch) "Scnp) To nb) "10 (np) "11(nb) "12 (nb) 


I 0.0430 0.6440 0.0112 0.2057 0.0070 Tey Tapp) Tm) TO ~My) 70-423 


3.1. Conditions for Convergence 

Convergence test is the method for testing the convergence of data collected prior to the previous 
one to get the roughly graph clearer. In Figure 6 shows the insignificant changing between the elements of 
two data even though it is known, that the Huffman code is the best sequence (13-length). Therefore, the 
multilevel sequence H(b) can be dignified as a good sequence. These sequences are: 

Huffman sequence: [1.0 1.081 0.597 0.79 0.686 -0.905 -1.086 0.905 0.686 -0.791 0.597 -1.08 1.0] 

Multilevel sequence [generated sequence H(b)]: [1.0000 0.2473 0.9478 0.4864 1.3280 0.3198 
0.0088 0.5767 -1.1750 0.4685 -0.9392 0.3623 -0.9817] 











o> Huffman Sequence 
=== Multilevel Sequence (good ACF) 

















oOo 








Sequences 






































6 8 
Length of Sequences 


Figure 6. Graph plotting for Huffman Sequences and Multilevel Sequences that are good in the ACF 


4. CONCLUSION 

This work has been presented methods for the design of multilevel sequences that have pulse like 
autocorrelation functions and relatively high energy ratio. This accomplishes to investigations into discrete 
coding techniques for improving range resolution and clutter performance of radar systems. The waveform 
considered in this work besides representing an interesting mathematical area, are also of practical 
significance in related fields such as sonar, navigation, and digital communications. 

The proposed methods contributes the generation of waveforms that have desirable properties for 
improving the range resolution of radar. Throughout this work digital processing has been assumed. The 
application of digital processing techniques to radar becomes more practical as compactness, cheapness and 
operational speed of digital micro circuits continue to increase. Although modern optical processing 
techniques sometimes provide an attractive alternative, the use of digital method with its inherent flexibility 
and reliability offers many advantages. To mention but a few, it simplifies pulse compression and real time 
multi-dimensional analysis of input data in range, Doppler, bearing, etc. Furthermore, it also offers 
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considerable advantages in post detection and display processing. In addition the use of digital processor will 
in many cases reduce future system modifications to easy and inexpensive software changes, rather than 
requiring costly hardware replacements. 

In conclusion, the multilevel code pulse trains design can be as the effective way to improve the 
overall performance of desired properties. It can enhance the radar performance by using the zeros pattern 
methods. It has the flexibility to dispose the noise ratio to improve the ER in that pulse trains. The ER valued 
is greater than original sequences but have poor characteristic of SLE compared between modified and 
without modification sequences but from the performance result it was efficient and flexible implementation 
compared the another methods, where it is the combination of ‘inverse filtering’ and ‘clipping’ provides a 
method of generating waveforms that have desirable properties for improving the range and Doppler 
resolution of radar [13-15]. 
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